Oscillations of Induced Magnetization in Superconductor-Ferromagnet 

Heterostructures 
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We study a change in the spin magnetization of a superconductor-ferromagnet (SF) heterostruc- 
ture, when temperature is lowered below the superconducting transition temperature. It is assumed 
that the SF interface is smooth on the atomic scale and the mean free path is not too short. Solving 
the Eilenberger equation we show that the spin magnetic moment induced in the superconductor is 
an oscillating sign-changing function of the product hd of the exchange field h and the thickness d of 
the ferromagnet. Therefore the total spin magnetic moment of the system in the superconducting 
state can be not only smaller (screening) but also greater (anti-screening) than that in the normal 
state, in contrast with the case of highly disordered (diffusive) systems, where only screening is 
possible. This surprising effect is due to peculiar periodic properties of localized Andreev states in 
the system. It is most pronounced in systems with ideal ballistic transport (no bulk disorder in the 
samples, smooth ideally transparent interface), however these ideal conditions are not crucial for the 
very existence of the effect. We show that oscillations exist (although suppressed) even for arbitrary 
low interface transparency and in the presence of bulk disorder, provided that /ir S> 1 (t - mean 
free path). At low interface transparency we solve the problem for arbitrary strength of disorder 
and obtain oscillating magnetization in ballistic regime (hr S> 1) and nonoscillating magnetization 
in diffusive one (/ir <C 1) as limiting cases of one formula. 

PACS numbers: 74.45.-l-c, 74.78.Fk, 74.78. Na 



INTRODUCTION 



Spin structures of the microscopic states of s-wave 
superconductors and ferromagnets are opposite to each 
other'^"^. Superconducting pairing interaction leads to 
formation of electron Cooper pairs with opposite projec- 
tions of spins, whereas the exchange field tends to align 
electron spins in the same direction. This counteraction 
on the microscopic level results in a competition between 
macroscopic superconducting and magnetic states. The 
suppression of the superconducting order parameter and 
the transition temperature by the exchange fieldS^ and 
the reduction of the magnetic spin susceptibility in the 
superconductors'''^ are well-known examples of this com- 
petition. 

Among suitable experimental systems for studying the 
interplay of the superconductivity and spin magnetism 
are superconductor-ferromagnet (SF) heterostructures 
(for reviews see Refs. @Q@)- Above the superconduct- 
ing transition temperature Tc, the superconductor is in 
its normal state, and the total magnetic moment Mtot of 
such system is given by the intrinsic magnetic moment 
of the ferromagnet MpQ. Below a magnetic moment 
M{h) induced by the presence of superconductivity ap- 
pears, and the total spin magnetic moment of SF system 
in the superconducting state is Mtot — Mpo + M[h). 
The induced magnetization may be caused by both the 
Meissner currents (orbital effect) and the spin polariza- 
tion (spin effect). If the sizes of the ferromagnet and the 
superconductor are small compared to the London pene- 
tration length, then the orbital effect is small compared 



to the spin effect^. In this work we assume this situation, 
since we want to study the effect related to the spin po- 
larization. Therefore M{h) is the induced spin magnetic 
moment throughout the paper. 

What direction of M{h) relative to MpQ one would 
expect? The above-mentioned competing behavior of su- 
perconducting and magnetic phenomena suggests that 
M{h) is opposite {M{h) < 0) to Mpa and thus reduces 
Mtot- In other words, the induced magnetization screens 
the intrinsic magnetization of the ferromagnet. The idea 
of spin screening of the ferromagnet 's magnetization by 
the superconductor in SF systems was first brought for- 
ward in Refs.iQj. In these publications the cases of a ferro- 
magnetic planar film and spherical grain were considered 
and it was shown that indeed M{h) < 0. 

Experiments carried out on various SF structures con- 
firm indirectly the idea of the screening proposed in 
Refs. H In Ref. [l3 a V-Pdi_:,Fe^ SF bilayered struc- 
ture was studied by a magnetic resonance technique and 
a 50%-decrease in Mtot was discovered as the temper- 
ature was lowered from T = Tc ~ 4K to T « 1.5K. 
In Ref. 1^ the neutron reflectometry was performed on 
multilayered SF structures consisting of ferromagnetic 
La2/3Cai/3Mn03 and superconducting YBa2Cu307 lay- 
ers. The obtained reflectometry spectra were discussed 
in context of the screening effect predicted in Ref. 0. 

In Refs. i2i the screening effect was studied under the 
following assumptions: i) "diffusive" limit (the mean free 
path of electrons I is much smaller than both the size 
of the ferromagnet d and the superconducting coherence 
length ^5); ii) the exchange field of the ferromagnet h 
is small compared to the Thouless energy i?Th = D/d?' 
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FIG. 1: (a) Oscillations of induced magnetic moment M(h) 
of a clean SF system as a function of h/eA ~ 2hd/vF (solid 
line shows —M{h)). The graph is plotted for the case of zero 
temperature T = and thin F- layer d <^ (^s- At hd <^ vp we 
have M{h) = —XNhd, indicating complete screening in the 
case itinerant ferromagnet. Dashed graph shows the function 
Z^{h/eA) (see Eq. Ill^l 'l. (b) The geometry of the system 
and spatial distribution of the density M{h, x) of the induced 
magnetization {h — hi - screening, h = h2 - anti-screening). 

(-D - difFusion coefRcient in the ferromagnet, throughout 
the paper we employ units, in which the Planck constant 
n = l): 

h^ETh, (1) 

and the ferromagnetic film is thin {d <ti ^s)- Condition i) 
allowed to use Usadel equation and condition ii) to treat 
the effect of ferromanget's exchange field /i as a pertu- 
bation in this equation. For the induced spin magnetic 
moment M{h) the following result was obtained: 

M{h)^~{XN~Xs{T))hd, (2) 

where Xs{T) is the magnetic susceptibility of a bulk su- 
perconductor, xn is the magnetic susceptibility of the 
superconductor in the normal state {xs(Tc) — Xn) (For 
exact expression for xs{T), see Eq. (|15l) '). The result 
Eq. ^ is quite universal^i, since it is independent of the 
strength of potential disorder and interface transparency. 

The following questions arise. First, how robust is this 
perturbative result Eq. (0) to the type of orbital electron 
dynamics, in particular, what result would be obtained 
in the opposite case of a clean ballistic system. Second 
and more interestingly, how does induced magnetization 
M{h) behave for sufficiently large exchange field, when 
its effect cannot be considered as a perturbation anymore 
and how the type of electron dynamics affects this behav- 
ior. The theory presented below shows that the behavior 
of induced magnetization M{h) in ballistic SF systems in 
nonperturbative regime can be very remarkable. 

First, let us define more precisely what we mean by 
"perturbative" and "nonperturbative" regimes for an SF 
system in general case. For a generic SF system with 
arbitrary bulk disorder in S and F regions and arbi- 
trary interface transparency an important energy scale 
is e* = 1/t* , where r* is the characteristic time spent by 



electron in the ferromagnet. (Ferromagnet is assumed to 
be of finite size d at least in one dimension). In ballistic 
system without or with relatively weak bulk disorder (the 
mean free path I ^ d) and with not too small interface 
transparency (t ^ 1) this energy scale is the Andreev 
energy: e* — — vp/d [vp is the Fermi velocity). In 
the case of low interface transparency t <^ 1 the time 
T* is enhanced due to the fact that electron has to hit 
the interface ^ 1/t times before it escapes from the fer- 
romagnet, therefore it stays in the ferromagnet ~ 1/t 
times longer compared to the case of good transparency. 
Thus, for the case of low interface transparency (t ^ 1) 
one gets e* = e^t for ballistic system. In the diffusive 
system (l <C d) with not too small interface transparency 
(t > l/d) e* = Ety, = D/d^ is the Thouless energy. 

Comparison of h with e* determines how strongly 
the exchange field h affects the spectrum of SF sys- 
tem compared to the spectrum of the corresponding 
SN(superconductor-normal metal) system with h = O.If 

h/e* < 1, (3) 

then the effect of exchange field is small and consequently 
physical quantities, such as the induced magnetization 
M{h), can be studied perturbatively. We refer to the 
regime as perturbative. On the contrary, when 

h>e\ (4) 

the spectrum of SF system is significantly altered by the 
presence of exchange field h, and we refer to the regime 
Q as nonperturbative. 

A very interesting property of Andreev spectrum of 
SF system in nonperturbative regime is its periodicity 
as a function of parameter h/e* with the period of the 
order of unity'^^ . (The physical reasons beyond this phe- 
nomenon and the origin of the conditions ©jO are given 
in Sec^l) This periodicity is known to reveal itself in 
the oscillations of the Josephson critical current Ic{h) in 
SFS junctionsi^ii^. One could expect to find these peri- 
odic features of Andreev spectrum in other macroscopic 
quantities, such as the induced magnetization M{h). In 
this paper we extended the analysis of Refs.Hjof induced 
magnetization in SF structures to nonpertubative case 
and found that in ballistic systems this is indeed the case. 

Namely, we considered an SF bilayered system with 
the ferromagnet of the thickness d and the superconduc- 
tor of the size much greater than its coherence length 
^5 (Fig. 1(b)). Solving the Eilenberger equation we find 
that in a system without or with relatively weak bulk 
{I ^ d) or surface disorder and with perfect SF interface 
transparency (t = 1) the induced spin magnetic moment 
M{h) is an oscillating sign-changing function (Fig. 1(a)) 
of parameter h/cA = 2hd/vp with the quasiperiod ap- 
proximately equal to tt. Therefore, the total magnetic 
moment 

Mtot = Mpo + M[h) (5) 

can be either smaller {M{h) < 0) or larger {M{h) > 0) 
than that in the normal state M^^^ = Mpo depending on 
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the value of the product hd. The oscillations of M{h) 
are most pronounced for the system with ideal trans- 
port properties: ballistic electron motion, perfect inter- 
face transparency. However, these ideal conditions are 
not crucial for the very existence of oscillations. To 
verify this we have considered the case of low SF in- 
terface transparency {t <^ 1) and arbitrary disorder in 
the ferromagnet, described by the scattering time t The 
limit of low SF interface transparency ^ 1) is use- 
ful from a methodological standpoint, since it allows to 
solve the Eilenberger equation for the system with arbi- 
trary strength of disorder. This allows one to study not 
only the limiting diffusive and ballistic cases, but also the 
crossover between them. It appears that the influence of 
disorder on the behavior of induced magnetization M{h) 
is governed by parameter hr. In the limit /it ^ 1 sign- 
changing oscillations of M (h) exist (the quasiperiod is 
still h* w t^£a), although their magnitude is suppressed in 
t as ^ and exponentially in d/l. On the contrary, in the 
opposite case hr ^ 1 we get that M{h) does not exhibit 
oscillations, being negative {M{h) < 0) for all h <C 1/t. 
The condition /ir <C 1 (together with I <^ d) corresponds 
to the "diffusive" limit of the Usadel equation and the re- 
sults obtained in this case from the Eilenberger equation 
can be recovered from the Usadel equation. 

We mention that nonoscillatory result for M{h) in the 
diffusive limit is in contrast with the behavior of the 
Josephson critical current Ic{h) in SFS junctions. Os- 
cillations of Ic{h) are not destroyed by disorder and per- 
sist (although exponentially suppressed in h/ETh) even 
in the "diffusive" limit hr ^ 1. The period of these os- 
cillations is h* ^ e* = Etyi. These oscillations were ob- 
served experimentally in Refs. Oj, 17 18 19-20., for further 
references see review articles Refs. Thus, oscillations 
of the induced magnetization M{h) turn out to be more 
sensitive to disorder than those of the Josephson critical 
current Ic{h). 

Our analysis shows that for moderate bulk disorder 
(I ^ d, I may also qualitatively include surface disorder 
of the interface) and not too small interface transparency 
(t ~ 1) the magnitude of oscillations is still quite no- 
ticeable, thus giving hope for experimental check of our 
predictions. Since for hr ^ 1 oscillations of M{h) are 
sustained, oscillatory behavior of M(h) should be attain- 
able even in the presence of disorder in the case of suffi- 
ciently strong ferromagnets. Since the exchange field h is 
hardly variable in the experiment, one may hope to ob- 
serve the oscillations of M(h) performing measurements 
on samples with different thickness d. We also note that 
the case of thin ferromagnetic films d <C is the most 
interesting for experiment: experimentally relevant ex- 
change fields are h Tc, one needs ~ e^i to observe 
oscillations, thus d/£,s ~ Tc/eA ~ Tc/h <C 1. 

As a limiting case of our analysis we obtain that in 
the clean case for small exchange fields {h <^ te^) and 
a thin ferromagnetic film {d <^ S,s) the induced magne- 
tization M{h) is given by the universal result Eq. ||2Jl. 
This complements the analysis of Refs. and suggests 
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FIG. 2: Semiclassical description of Andreev spectrum of an 
SF system. Each classical trajectory that starts and termi- 
nates at SF interface corresponds to a set of discrete Andreev 
energy levels, which are obtained from the Bohr-Sommerfeld 
rule Eq. 

that this result holds in perturbative regime in SF 
systems with arbitrary strength of potential disorder. 

The paper is organized as follows. In Sec.^we provide 
qualitative quasiclassical description of Andreev spec- 
trum in SF systems and show how the conditions 
and the periodicity of the spectrum arise. In Sccs. lIIll we 
consider the limit of clean sample and perfect interface 
transparency, when the oscillatory behavior of induced 
magnetization is most pronounced. We present the sys- 
tem and formalism of the Eilenberger equation used to 
derive the expression for the induced magnetization and 
analyse this expression in detail. The connection between 
the predicted effect and the properties of Andreev spec- 
trum is discussed. In Sec. IIVI we show that our assump- 
tions about ideal transport properties (perfect interface 
transparency, ballistic electron motion) of the system are 
not crucial for the existence of oscillations. We present 
the results for the case of low SF interface transparency 
and disordered ferromagnet and show that oscillations 
exist in such limit provided some conditions on parame- 
ters are met. Finally, we conclude with Sec. In Ap- 
pendix the general formulas for the case of arbitrary in- 
terface transparency and clean samples are given. 



II. QUALITATIVE PHYSICS OF ANDREEV 
SPECTRUM IN SF SYSTEMS 

The spectrum of a superconductor-ferromagnet(SF) 
(or superconductor-normal metal(SN)) system is given 
by Andreev states, which are the states of electron-hole 
pairs localized in the F(N) region due to Andreev reflec- 
tion at SF(SN) interface. A qualitative understanding of 
the properties of Andreev spectrum in SF systems can be 
obtained from the following semiclassical picture (Fig.l^J. 
We assume the case of ideal SF interface transparency 
here for simplicity. (Qualitative analysis of the present 
section generalizes the discussion of Andreev spectrum of 
SN systems done in Sec. II of Ref. |^ to SF systems.) 

First, consider an SN system without exchange field. 
Suppose an electron with an energy e relative to the 
Fermi level tp inside the superconducting energy gap A 
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(|e| < A) is travelling along some classical path in N re- 
gion. If the electron hits SN interface, it is reflected as 
a hole with the energy — e. A peculiar property of An- 
dreev reflection is that the momentum of reflected hole is 
opposite (apart from a small angular mismatch ~ e/ef) 
to that of the incident electron. Therefore the reflected 
hole will travel along the same path as the incident elec- 
tron but in the opposite direction. If this path hits the 
interface again, the hole is reflected back as the elec- 
tron. Within the quasiclassical Bohr-Sommerfeld(BS) 
approach if the total action of such process is an inte- 
ger multiple of 27r?i, then such electron-hole pairs forms 
a bound state. 

Thus, each classical path that starts and terminates 
at SN interface corresponds to a set of discrete Andreev 
levels, energies of which can be obtained from the BS 
rule. For each such path 7 of the length the action of 
the electron traversing this path is S{e) — p{e)L^, where 

p{e) = ^2m(e_F + e) ^ pp + e/fF is the absolute value 
of electron's momentum. The action takes the form 

S{e) = pfL^ + er^ 

where = L^/vp is the time of traversal of path 7. The 
action of the hole (a missing electron) is — S'(— e) and the 
contributions from Fermi length scales cancel each other. 
The BS rule gives 

S{e) - S{-e) - 20(e) = 2er^ - 2(/)(e) = 27rn, 

where 0(e) = arccos(e/A) is the phase of Andreev reflec- 
tion and n is integer, and thus discrete Andreev levels 
e = e^(n) corresponding to path 7 are determined from 
the equation: 

er^ = nn + <p{e). (6) 

Generalization to SF systems is straightforward 
(Fig. 121) . In SF system one should distinguish between 
Andreev states with electron spins directed along (t) and 
opposite to (J,) the exchange field h. "Up" states e| ac- 
quire a shift —h and "down" states e| acquire a shift +h 
in the F region. The BS rule reads(-|- and — correspond 
to t and I, respectively): 

5(eT,x ±h)- 5(-(et,x ± h)) ~ 20(et,x) = 27rn 

Note that the argument of (f>{ei^i) does not acquire the 
shift ±h due to exchange field, since exchange field is ab- 
sent in the superconductor and therefore does not affect 
the process of Andreev reflection. Thus, the equation for 
Andreev levels in SF system reads: 

(e|,j, ± /i)r^ = 7rn + 0(e|,|) (7) 

The time r* defined in Sec. |l]as a time spent by elec- 
tron in the F region is the characteristic time for t^. If 
hr* <C 1, then Andreev levels of SF system obtained from 
Eq. O are only slightly different from those of the cor- 
responding SN system with h — obtained from Eq. 



This condition corresponds to perturbative regime 
If hr* > 1, then Eq. Q is significantly different from 
Eq. ® and this condition corresponds to nonperturba- 
tive regime Further, one sees that Eq. {Tj) is invariant 
under the periodic translation (/ir^ —^ hTy + n). Thus the 
solutions of Eq. {Tj) are periodic: 

This periodicity of Andreev spectrum should reveal itself 
in the behavior of macroscopic quantities of SF system 
as a function of parameter hr* . This is indeed the case 
for Josephson critical current in SFS junctions and, as 
we show in this paper, for the induced magnetization. 

III. BALLISTIC CASE, IDEAL SF INTERFACE 
TRANSPARENCY 

A. System and method 

We start our analysis with the case of ideal ballistic 
electron transport: absence of bulk disorder and perfect 
interface transparency t = 1. In this case the oscillations 
of induced magnetization are most pronounced. 

We consider an SF bilayered system pictured in Fig. 
1(b): the superconductor occupies the half space x > d 
and the ferromagnetic layer (F-layer) of the thickness d 
is located in the region < x < d . We assume that 
there is no bulk disorder neither in the F-layer nor in the 
superconductor and that the interfaces are ideal, namely, 
the superconductor- ferromagnet(SF) interface {x = d) is 
perfectly transparent and the ferromagnet-vacuum (FV) 
interface {x = 0) is specular. 

We study the problem solving the Eilenberger 
equations^ for the quasiclassical Green's function 
g{oj,n,r). Here oj — 7rT(2m -I- 1) is the fermionic Mat- 
subara frequency (T is the temperature and m is an in- 
teger), n is the unit vector representing the direction 
of the electron momentum on the Fermi surface and r 
is the radius vector. Due to the geometry considered 
g{u!,n,r) = g{Lu,nx, x), where is the projection of n 
on the a;-direction, —1 < < 1- The Green's function 
g{uj,nx, x) is a matrix in the tensor product of Gor'kov- 
Nambu and spin spaces: 

g ^ gi°Ti+ 52 0T2+ .93 o T3, 



g^^g°i + gta,. (8) 

Here r^, i — 1,2,3 are the Pauli matrices in the Gorkov- 
Nambu space, and 1, ct^ are the unity and Pauli matrices 
in the spin space {z denotes the direction of the exchange 
field in the spin space). Diagonal representation ^ in 
the spin space is possible in the case of homogeneous 
magnetization, which is assumed here. 

The quasiclassical approach implies that the following 
conditions are satisfied: h <C ep, ppd 3> 1, where ep, pp 
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are the Fermi energy and momentum. However, we stress 
that these conditions are required only for the appUca- 
biHty of the method used, but not for the very existence 
of the oscillations. The results are still valid qualita- 
tively in the case of strong ferromagnet of atomic thick- 
ness {h < ep, ppd 1), although the effect is reduced. 
The same concerns possible mismatch of electronic prop- 
erties (such as density of states vp, Fermi velocity vp) in 
the ferromagnet and superconductor. It is assumed here 
that they are the same. 

The density of the induced spin magnetization can be 
expressed in terms of the quasiclassical Green's function 
in the following way: 

M{h,x) = A^s((^f(x)^T(^)>-(V'i"(^)^i(^)) = 
= —fiB>^pTy2 I dri^ g^{uj,n^,x), 

where /xb is the Bohr magneton and vp is the electron 
density of states at Fermi level per single spin projection. 
The total induced magnetic moment of the system (per 
unit square in the plane parallel to SF interface) is ob- 
tained by the integration of M{h, x) with respect to x: 

r+oo 

M{h) = / dxM{h,x). (9) 
Jo 

The relation between the intrinsic magnetic moment of 
the ferromagnet MpQ and the exchange field h acting on 
the electrons depends on the model of the ferromagneU^. 
Usually MpQ = Mp^^i + Mpoioc is combined of the con- 
tribution Mpoei = XN^d produced by free electrons and 
the contribution MpQioc = axN^d produced by localized 
magnetic moments (a is some phenomenological constant 
and xn is the normal metal spin susceptibility). In the 
case of itinerant ferromagnet the localized magnetic mo- 
ments are absent (a = 0) and Mpo = Mpo^i = XNhd. 

We emphasize that magnetic moment M{h) Eq. @ 
expressed in terms of the quasiclassical Green's function 
is the magnetic moment induced by the presence of super- 
conductivity in the system. It is determined by the prop- 
erties of the energy spectrum on the scale ~ Tc near the 
Fermi level. It does not include the part MpQ^i = XNhd 
of the intrinsic magnetic moment Mp^ of the ferromagnet 
produced by free electrons that originates from the en- 
ergy shift of the entire electron band. The latter cannot 
be taken into account within the quasiclassical approach 
and should be added separately. The total magnetic mo- 
ment of the system is given by Eq.®. 



The Eilenberger equation for the system without dis- 
order reads: 

vpn^d^g + [{uji - ih{x)a^) o T3 + A{x)i o T2,g] = 0, 

(10) 

where h(x) is the exchange field in energy units, A(a;) 
is the superconducting order parameter, and brackets [, ] 
stand for the commutator. The exchange field is con- 
tained in the F-layer only and assumed to be constant 
within the layer: h(x) = h, if < x < d, and h{x) = 0, 
a X > d. We assume there is no BCS interaction be- 
tween electrons in the F-layer and therefore we always 
have A(a;) = for < a; < d. 

In principle, the exact order parameter A(x) has 
to be found self-consistently from the solution of the 
Eilenberger equation (|10|l supplemented by the self- 
consistency equation for the superconducting order pa- 
rameter. The order parameter A(a;) approaches a bulk 
BCS value A at large distances from the F-layer x ^ £,s, 
but is partially suppressed near the F-layer. Computa- 
tion of A(a::) self-consistently is a hard analytical prob- 
lem. Fortunately, our main result about the oscillatory 
sign-changing behavior of AI (h) is not sensitive to the ex- 
act shape of A(a;). Therefore, we assume that A{x) = A 
for all X > d and perform calculations under this assump- 
tion. 

Eq. H1U|I must be supplied by proper boundary condi- 
tions at SF and FV interfaces'^. The limit of the ideal 
transparency of SF {x ~ d) interface implies that the 
Green's function is continuous: 

g{u!, n,j;,x = d-0) ^ g{uj, n^.x = d + Q). 

At FV {x — 0) interface the specular reflection condition 
reads: 

g{uj, n^;, X = 0) = g[uj, -n^, x = 0). 
At X ^ ^5 the solution approaches the BCS bulk result: 

z Q CiJ 

9t = 5l = 0, g^^ fs ^ —===, g^^ gs ^ —===. 

V + A^ VtjJ + A^ 



B. Analysis 

Under the made assumptions the solution of Eq. H10() 
is straightforward and we obtain: 



g3{uj,n^,x) 



o-IS 



sm2H J 1, < x < d, 

[cosh 17 + \gs\ sinhflf cos' H + [sinhSl + \gs\ cosh 17]^ sin' H \ cxp {-^^^^^^{^ 



X > d. 
(11) 
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Here n = \uj\/{eA\na:\), H = h/{eA\nx\) and ea = 
vf/ {'2d) is the Andreev energy: Ti/eA is the time the elec- 
tron travels from SF interface and back within F-layer 
with the velocity perpendicular to the interface. 

Inserting Eq. into Eq. © and integrating over x 
one obtains M{h). The key point of our analysis is that 
ig^{u!,nx,x) is a periodic sign-changing (note sin2iJ in 
the numerator in Eq. function oi H ~ h/{eA\nx\) 

and depends on the exchange field solely via this param- 
eter. 

General properties of M{h) can be summarized as fol- 
lows: 1) M{h) depends on the strength of exchange field 
solely via the combination h/eA — 2hd/vF', 2) for any 
temperature T < Tc and any ratio rf/^s the induced 
magnetic moment M{h) is an oscillating sign-changing 
function of h/eA with a quasiperiod h* ftA ~ tt, the 
amplitude of the oscillations decays monotonically as 
hjeA increases; 3) hence, M{h) can either have the same 
(M(/i) > 0) or opposite (M(/i) < 0) direction as Mpo, 
depending on h/eA', 4) at h/eA I we get AI{h) < 0, 
which indicates the screening of Mpo] 5) the magnitude 
of oscillations of M{h) is largest at T = and decreases 
as T increases; M{h) = at T = Tc. 

The spatial dependence of M(/i, x) shown in Fig. 1(b) 
is governed by g^{uj,nx, x): M{h,x) is constant within 
the F-layer and decays exponentially over the distance 
into the superconductor. If d ~ ^s, then parts of 
M(h) located in the F-layer 



dxM{h,x) 



and in the superconductor 



Msih) = 



+ 00 



dxM(h,x) 



are of the same order. If d 3> ^5, then the induced mag- 
netic moment is located predominantly in F-layer and 
M{h) w Mp, Ms/Mp ^ is/d < 1. In the opposite limit 
d <^ the induced magnetization is located mainly in 
the region of the superconductor of the size near the 
F-layer and M{h) « Ms, Mp/Ms - d/^s < 1- 

Below we concentrate on the experimentally more rel- 
evant situation of a thin F-layer d <C and illus- 
trate the announced properties of M{h) explicitly for 
this particular case. In this regime, the expression for 
M{h) « Ms{h) can be reduced to the form: 



M{h) 



-d^Bl^F^AT^Ty^^ 



A2 



Vt^2 + A2 

sin2i/ 



;2-H A2 cos^iJ' 



(12) 



First, we consider the zero temperature limit T — Q. 
Replacing the sum over lj by the integral T . . . ~ 
duj/ (27r) . . ., we obtain: 



M{h) = -2diiB>^peA I -^Z^{—t 



(13) 



where Zt^^x) = x, if — 7r/2 < x < 7r/2, and periodically 
continued to all x (linear "zig-zag" -type function with a 
period tt). The integral with respect to nx can easily be 
calculated and we obtain the function shown in Fig. 1(a). 

Close to superconducting transition point {{Tc — 
T)/Tc < 1) we obtain 

M{h) = -dfiB'^peATTTY^ — J — sin 

We see that M{h) is again an oscillating sign-changing 
function of h/eA with quasi-period h* /eA ~ ti", although 
the amplitude of oscillations is parametrically smaller 
than that at T = by (A(r)/Tc)2 - (T^ - T)/Tc < 1. 



C. Perturbative regime (h <^ teA) 



In the limit h <^ eA from Eq. (|12|l we get: 

(tj2 + A2)3/2 



M{h) = -xnttTY, 



hd (14) 



where xn — 2^bvp is the bulk spin susceptibility of the 
normal metal. Since 



Xs{T) = 1 - '^T^ 



A2 



{uj-^ + A2)3/2 



(15) 



is the bulk spin susceptibility of the superconductor, 
Eq. H14|) can be rewritten in the form Eq. From 
the formulas given in Appendix one obtains the same 
result in the case of clean samples and low SF inter- 
face transparency {t <C 1), provided that h <C teA- 
The total magnetic moment produced by free electrons 
is Mpoei + M{h) = xs{T)hd. At zero temperature 
Xs{T = 0) =0 and M{h) = -Mpoei = XNhd, i.e. the 
induced magnetic moment M{h) totally screens the part 
Mpoei of the intrinsic moment Mpo produced by free 
electrons. And interesting feature is that Mpoei is lo- 
cated in the ferromagnet, whereas M{h) is spread over 
the distance ^5 from the F-layer in the superconductor. 
Since for h eA and d <^ the exact order parameter 
A (a;) is only shghtly suppressed due to the ferromag- 
netic proximity effect, this result is justified even if the 
self-consistency condition for A(a;) is taken into account. 

As the same result Eq. was obtained in the oppo- 
site diffusive limit for h <C -Exh, we make a conjecture 
that Eq. (|2J) holds for arbitrary strength of potential dis- 
order, provided the general condition Eq. Q is met. The 
universality of result Eq. Q is reminiscent of the prop- 
erties of the bulk linear spin susceptibility of the super- 
conductor xs{T) (Eq. (|15|) '). It is also independent of the 
strength of potential disorder—. 



D. Andreev states 

The oscillations of induced magnetization are closely 
related to the properties of the energy spectrum of lo- 
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calized Andreev states in the systemic. The equation 
for Andreev energy levels e|.| with the electron's spin 
having the same (t) and opposite (|) direction as the 
exchange field (corresponding to + and — signs, respec- 
tively) reads: 



eta ±^ 



arccos 



A 



7rn, 



(16) 



where n is integer. Note that in the case of ballistic 
system this equation exactly coincides with Eq. used 
for qualitative considerations. Eq. is invariant under 
the periodic translation h/{eA\nx\) h/ {eA\nx\) + Trk {k 
is integer) in the same fashion as the Green's function 
is periodic in H. In the limit c? <C there exists 
only one level for a given rix and projection of spin2£: 



e|,4(iJ) = ±Acosi/, i/e[0,7r], 



(17) 



and periodically continued to all H. The states with 
H e [0,7r/2] + nn (e| > 0, < 0) contribute to the 
screening of Mpo (sin2iJ > 0, see Eq. (|12ll ). whereas the 
states with H e [7r/2,7r] + 7m (e| < 0, > 0) give rise 
to the anti-screening of Mpo (sin 2H < 0) . Due to the 
property 

e|,x(7r/2 + nn + 5H) = e|,|(7r/2 + mi - 6H), 

6H e [0,7r/2], such "up" and "down" states interchange 
in the energy space but since the spin direction is "at- 
tached" to them explicitly, this results in the opposite 
signs of contributions to M(/i). 



E. Self-consistency of order parameter A(a::) 

As it has been mentioned, the oscillating behavior of 
M{h) is insensitive to the exact shape of A(x) and there- 
fore persists if the self-consistency of A(a::) is taken into 
account. This is the case, because the periodic functions 
of H in Eq. (|ll|l arise from the solution of Eq. H10|) in 
the F-layer, where A(a;) — and the general solution 
can always be found explicitly. Due to this fact ig^ is a 
periodic sign-changing function of H and, hence, M{h) 
is a quasiperiodic sign-changing function of h/eA, inde- 
pendently of the exact shape of A(a;). Since for arbitrary 
T the proof is cumbersome, we illustrate it here in the 
simplest case, when T is close to Tf,. In this limit one can 
obtain the solution to Eq. (|10|l as an expansion in A (a;). 
In the lowest (quadratic) order we find for x > d: 

i f h 1 \ 1 

g^{u),nx,X) = --rrrrSm I 2 ■; r I — X 



2e V K\j 

A{x')e'T dx' A{y)e~tdy, 



where ^ = ^2\2\^ ■ This yields the form 

M{h) - ^ dt sin (^2^t^ F{t, h), 



where F{t^ h) is a positive monotonically decreasing with 
respect to t envelope function. The above integral can 
be both positive and negative depending on the value of 
h/eA- This is especially clear for h/eA ^ 1 when one can 
integrate by parts to obtain 



M{h) 



cos {2h/eA) 
2h/eA 



F{l,h). 



Therefore the induced magnetization is an oscillating 
function of the parameter 2h/eA regardless of the ex- 
act form of A(x). 



IV. LOW SF INTERFACE TRANSPARENCY 

Our assumptions about ideal transport properties of 
the system (perfect SF interface transparency, ballistic 
electron motion in the samples) are not crucial for the 
existence of oscillations of induced magnetization. The 
oscillations of M{h) exist for arbitrarily low SF-interface 
transparency (see Appendix). Moreover, they can exist 
in the presence of bulk disorder in the sample. 

To illustrate that we turn to the case of low SF inter- 
face transparency t = tijix) <C 1 (tijix) is a transmit- 
tance coefficient, see Appendix). In this limit the prox- 
imity effect is weak and one can take the effect of disorder 
into account by linearizing the Eilenberger equation with 
collision term^"*-^^. We assumed that superconductor is 
clean and ferromagnet is disordered and described by the 
mean free path I and scattering time t = l/vp- 

The Eilenberger equation in the F region {0 < x < d) 
for the "up" component g] = g reads (we omit the index 
"f here for brevity): 



vpnxdxg 



{uj-ih)T3 + —{g),g 
It 



= 



(18) 



where (5) = 1/2 Jj'^j^ dnxg{nx, x) is the angular averaging. 
In the S region {x > d): 



vprixOxg + [uJT3 -\- At2, g] = 0. 



(19) 



In the zcroth order in interface transparency {t{nx) = 0) 
the superconductor and ferromagnet are not linked and 
the solution is: 

g^p^ = sgnwTa, < x < d. 



-(0) 

9 s 



fsT2 +gsT3, X > d 



Next, we present the Green's function in the ferromagnet 
in the form 

9 = 9f^ + ^9, Sg = SgiTi + Sg2T2 + SgsTs, 

and, leaving only linear in 5g terms in Eq. (|18|l . arrive at 
the following equation for Sg: 



vpnxdxSg- 



., , sgnw 



sgn w 
2t 



[{Sg),T3]^0. 
(20) 
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First we need to obtain a linear in t solution for dg2 in 
the F region. From Eq. (|18|l we get: 



(21) 



where a^^; = 1 + 2(|a;| — i/isgna;)T. Boundary condition 
with vacuum reads: 



d.Jg2{x^0)^0 



(22) 



At the SF interface one must use Zaitsev boundary con- 
ditions IjAip for nonideal interface transparency. In the 
limit t ^ I one can expand them in t. First order in t 
gives: 

5gi{n.j„x = d - 0) = i^sgnujsgnn^fs. 
Using the T2-component of Eq. (|20|l 

lnxdxSg2 + ioiuj sgnuSgi — 0, 

we arrive at the boundary condition for Sg2 at SF inter- 
face: 



dx5g2{x = d - 0) = ^^y^/sac 



(23) 



Eq. H21|) must be solved for x G [0, d] with boundary con- 
ditions Eq. lf^ . l|^ . Condition Eq. allows to sym- 
metrically continue 5g2 to [— c?, 0] interval. We perform 
the Fourier transformationSiS^: 



-l-oo 



g2ix)^ J2 Sg2in)e^'-\ 



Next we solve Eq. (|19|l in the S region and get for 

X- - -(0) 
^9 = 9-93 ■ 



Sg ^ c{nx){sgnnxTi- gsT2 + fsT3)e "-fI"-! (26) 
where c{nx) is a symmetric (yet unknown) function of 

It follows from Eqs. (P|). 12t)|) that the induced magnetic 
moment is given by (again we assume the case of a thin 
F film {d <^ £,s) and therefore the induced magnetic mo- 
ment M{h) « Mg{h) is located in the superconductor): 

M{h) = ^BVp'^T^Ty^ I dux -—-^ =fslrnc{nx) 
^ Jo 

Expanding Zaitsev's boundary conditions to the second 
order in t, we obtain: 



Imc(n3: 



and therefore 



t{nx) 



gsluiSg2{nx,x = d ~ 0) 



fs9s 



M{h) = -fiBi^FVpT^Ty^^ -y= 

X / dnxnxt{nx)lixi6g2{x = d). 
Jo 

The needed quantity is 

oo oo 

6g2{x^d)^ ^ (-l)"<5g2(n) = ^ Sg*2{r 



n— — QO 



n— — oo 



where A:„ = im/d {n is integer) and Fourier coefRcients 
are 

1 



('^) = TTl / 5g2{x)e-^''"^dx. 



2d 

Calculating Fourier coefhcient of both sides of Eq. H21|l 
and taking the boundary condition H23|) into account, we 
get: 

592{n) = ^ m2{n)) + {-lTeAT\nx\t{nx)fs) , (24) 

where L„ = Pn^.k'^ + a^. Angular averaging of Eq. 
gives: 



{592{n)) 



(-l)"«./5eAr(^) 



and thus 



5g2{n)^{~ir5gl(n), 



' |"x|t\ 



592 (n) - ^e^rfs + M \ . (25) 



Below we analyse two different limiting cases depend- 
ing on the strength of disorder. It appears that the in- 
fluence of disorder on behavior of induced magnetization 
M{h) is governed by parameter hr rather than d/l. 



A. Quasiballistic case (Jit 2> 1) 

If hT ^ 1, then a^^/L^ ^ ^jci^ ^ VC*''') 1 ^.nd one 
can neglect the first term in parentheses in Eq. (|25|l and 
get: 



Summing the series, we get 

. , t / \uj\ ~ ih sgn uj + 
Sg2{nx,n) = -/scoth ■ — ■ ^ 

and 



2 _L ^2 

t^l — ihsgnu) -t- ^ 



X / duxUxt (ria;)Im coth 

Jo V eA|n; 
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For not too small disorder {d/l > 1) (and a thin F-layer 
{d < (,s) we get: 



A2 



M{h) = -\t^B'yFVFTrTY, , 2 ^ ^2)2 

X / drLxUxt^ {rix) eyi]> ( ) sin2i7. 

Jo V "^xl J 

One sees that oscillatory behavior is sustained, although 
the magnitude of oscillations is suppressed in t{nx) <C 1 
and 1 1 d. Extrapolation of this formula to not too small 
i ^ 1 gives that for moderate disorder / > d the magni- 
tude of oscillations is quite comparable to the ideal case 
Eq. 1)12(1 . We also mention here, that spin-orbit scatter- 
ing should also suppress oscillations of M{h). Spin-orbit 
scattering is neglegible, if the corresponding mean free 
path Iso 3> d. This condition is always satisfied for / ^ d, 
because Iso^ I- 



B. Diffusive case {hr <C 1) 

If /ir <C 1 and I ^ d, then a^{l/Ln) — + 1 and the 
main contribution to Sg2 (n) comes from the first term in 
parentheses in Eq. H25|l . which has a (diffusion) pole, and 
the second term can be neglected: 



<5.92(") = ^Afs{\nx\t) 



1 



2{\uj\ - ihsgnuj) + Dkl 



where D — vpl/i is the diffusion coefficient. Summing 
the series, we get 



5g2{x = d) = fs{\nx\t) 



cothA/ '^l^l^^^fg""^ 



2/ 



2[\iLj\~ih sgn uj) 



and 



M{h) — — ^B^pvpT^T 



fs9s 

^ VA2 -I- w2 



coth 



2{\u\ — ih sgn u;) 



2{\Lo\ — ih sgn oj) 



where £^Th = i'/d^ is the Thouless energy. This result 
is valid, if < <C Ijd and h ^ ie^. Since the conditions 
/it <C 1 and / <C d correspond to "diffusive" regime, this 
result could be obtained from the Usadel equation. 

Interestingly, Im bg2 (x = d) does not oscillate as a 
function of h, even though it contains trigonometric func- 
tions. Therefore, we obtain that in the diffusive limit 
/it <C 1 the induced magnetization M{h) is not oscilla- 
tory and always negative. 



V. CONCLUSION 

In conclusion, we have shown that in SF systems the 
total spin magnetic moment in the superconducting state 
can be both smaller and larger than that in the normal 
state. The effect is due to peculiar periodic properties of 
Andreev states in SF systems that result in oscillatory 
sign-changing behavior of the superconductivity-induced 
magnetization of the system. The predicted effect is ex- 
pected to be best observable in relatively clean SF sys- 
tems with good quality of interfaces. Practically this 
means that the mean free path I should be larger than 
the "exchange length" /gxc = vp/h. This condition can 
be fulfilled in the case of sufficiently strong fcrromag- 
nets. On the other hand / should not be much smaller 
than the thickness of the ferromagnetic film d. We ig- 
nored a change in the magnetic moment M{h) caused by 
the Meissner currents assuming that the thicknesses of 
the ferromagnet and superconductor are smaller than the 
London penetration length. In this case the contribution 
of these currents to M{h) is small. Spontaneous orbital 
effects in clean SF structures were studied in Ref. 
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erostrukturen and German- Israeli Foundation for Scien- 
tific Research and Development for financial support. 



APPENDIX A: CLEAN SAMPLES, ARBITRARY 
SF INTERFACE TRANSPARENCY 



Boundary conditions for nonideal interface trans- 
parency have been derived by Zaitse\4^. They are ex- 
pressed in terms of the antisymmetric 

d{nx,x) = {g{nx,x) - g{-nx,x))/2 
and symmetric 

s{nx,x) = {g{nx,x) -\- g{-nx,x))/2 

parts of Green's function in the following way: 

d{{l - t){s+ + s^f + {s+- s^f) = /(s+-s_.)(s+-Ks_). 

(Al) 

Here d — d{nx,d), s+ = s{nx,d 0), s_ — s{nx, d — 0). 
Antisymmetric part a{nx^ d) is continuous at the bound- 
ary X = d. The transmittance coefficient tijix) can vary 
from t{nx) — for nontransparent interface (e.g. bound- 
ary with vacuum) to t{nx) = 1 for perfectly transparent 
interface. 

In the case of clean samples and arbitrary transparency 
t{nx) one must solve Eilenberger equation (|10|l with 
boundary conditions l|Aip . We obtain (a; > 0): 
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sinh(f7 - iH) 



[gs cosh(17 - iH) + {2/t{n^) - 1) sinh(f7 - iH)f + /| 



exp ^ — — A2 



in the superconductor {x > d) and 



5S cosh(ri - iH) + (2/i(n^) - 1) sinh(f7 - ii/) 
[5s cosh(17 - iH) + {2/t{nx) - 1) sinh(f7 - ii?)]^ + /| 

I 



(A3) 



in the fcrromagnet (0 < x < c?) (For notation, see 
Secs. lIII Aim B|) . One can check that induced magnetiza- 
tion M (h) foUowing from these formulas is an osciUating 
sign-changing function for arbitrary ^(n^:), < t{nx) < 1- 



The magnitude of osciUations of M{h) is greatest at 
t{nx) — 1 and decreases as t{nx) decreases; M{h) = 
at t{nx) = 0. 
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